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Abstract — Explicit Finite-Difference method was used to obtain the solution of the system of the non-linear ordinary 
differential equations which transform from the non-linear partial differential equations. These equations describe the steady 
magneto-hydrodynamic flow of an oldroyd 6-constant non-Newtonian nano-fluid through a non-Darcy porous medium with 
heat and mass transfer. The numerical formula of the velocity, the temperature, the concentration, and the nanoparticles 
concentration distributions of the problem were illustrated graphically. The effect of Darcy number D a , Forchheimer number 
F s , magnetic field parameter M, local temperature Grashof number G n local nanoparticle Grashof B n Prandtl number P n 
Dufour number N d , Brownian motion parameter N h , Thermophoresis parameter N„ Lewis number L e , Sort number L d , 
Chemical reaction parameter R c , and Chemical reaction order m on those formula were discussed at the values of material 
parameters (a ± = 8 , a 2 = 2) specially in the case of pure Coutte flow. Then, the effects of modified pressure gradients on 
those formulas were discussed in the case of pure Poiseuille flow and the generalized Couette flow. Also, an estimation of the 
global error for the numerical values of the solutions is calculated by using Zadunaisky technique. 

Keywords — Finite-difference method, Heat and Mass transfer, MHD flows, Non-Darcy Porous medium, Oldoryd di- 
Constant non-Newtonian nano-Fluid. 


I. Nomenclature 


A 

Chemical Reaction rate constant 

Be 

Chemical reaction parameter, defined by Eq. (42) 

B, 

Local nanoparticle Grashof number. 

S 

The dimensionless concentration 

c 

defined by Eq. (42) 

Forchheimer number 

t 

The time 

C 

The concentration of the fluid 

T 

The fluid temperature 

Cl 

The concentration at lower plate (y = 0) 

7) 

The temperature at lower plate (y = 0) 

C 2 

The concentration at upper plate (y = //) 

t 2 

The temperature at upper plate (y = II) 

C s 

Nanoparticle susceptibility 

V 

The velocity vector 

D a 

D 

B 

D 

T 

Darcy number, defined by Eq. (42) 

Brownian diffusion coefficient 

Thermophoretic diffusion coefficient 


Greek symbols 

E 

Electrical field 

<p 

The nanoparticles phenomena 

F 

The external force 

<t> 

The dissipation function 

F s 

Forchheimer number, defined by Eq. (42) 

0 

The dimensionless nanoparticles 

G 

Gravitational acceleration 

9 

The dimensionless temperature 

G r 

Local temperature Grashof number, 
defined by Eq. (42) 

CC i f CX-2 

The non-Newtonian parameters , defined by Eq. (36) 

H 

J 

The magnetic field 

The current density 

V 

Gradient operator 

k 

Thermal conductivity 

V 2 

Laplacian operator 
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k * 

Permeability constant 


Kj 

Thermal diffusion ratio 

Pe 

L e 

Lewis number, defined by Eq. (42) 

V 

La 

Sort number, defined by Eq. (42) 

c v 

M 

Chemical Reaction order 

pf 

M 

Magnetic field parameter, defined by Eq. 
(42) 

Pp 

N h 

Brownian motion parameter, defined by Eq. 
(42) 

(/r)/ 

N d 

Dufour numbe, defined by Eq. (42) 

(pc)p 

N, 

The thermophoresis parameter, defined by 


Eq. (42) 

(7 

P 

The fluid pressure 

T 

P, 

Prandtl number, defined by Eq. (42) 

CCj 1 t CCq 


the dynamic viscosity of fluid 
the magnetic permeability 
The kinematic viscosity (”/p ) 

The specific heat capacity at constant pressure 
The density of the fluid 

The density of the particle 
heat capacity of the fluid 

effective heat capacity of the nanoparticle material 

Electrical conductivity of the fluid 
The Cauchy Stress tensor 

Volumetric thermal and solute expansion coefficients of the 
base fluid 


II. Introduction 

The study of non-Newtonian fluids has gained much attention recently in view of its promising applications in engineering 
and industry. Such fluids exhibit a non-linear relationship between the stresses and the rate of strain. Due to non-linear 
dependence, the analysis of the behavior of the non-Newtonian fluids tends to be much more complicated and subtle in 
comparison with that of Newtonian fluids. Flow of fluids with complex microstructure (e. g. molten polymer, polymer 
solutions, blood, paints, greases, oils, ketchup, etc.) cannot be described by a single model of non-Newtonian fluids. Many 
models that exist are based either on natural modifications of established macroscopic theories or molecular considerations. 
In general, the equations of motion for non-Newtonian fluids are of higher order than the Navier-Stokes equations and thus 
one need conditions in addition to the usual adherence boundary condition [28], 

Guillope and Saut [2] has established existence results for some shearing motions of viscoelastic fluids of Oldroyd type. 
Some exact solutions of an Oldroyd 3-constant fluid are studied in [4], [5], [11], [13]. Baris [21] investigated the steady flow 
of an Oldroyd 6-constant fluid between intersecting planes using the series expansion method. Hayat et al. [23] studied the 
Couette and Poiseuille flows of an Oldroyd 6-constant fluid with magnetic field by using the Homotopy analysis method. 
Hayat et al. [25] studied the steady flow of a magneto -hydrodynamic Oldroyd 6-constant fluid by the motion of an infinite 
plate using the Homotopy analysis method. Wang et al. [30] investigated the non-linear magnetohydrodynamic problems of 
an Oldroyd 6-constant fluid by using analytical method and the finite-difference discretization method. Hayat et al. [24] 
studied the effect of the slip condition on flows of an Oldroyd 6 -constant fluid. Rana et al. [14] studied the Hall effects on 
hydromagnetic flow of an Oldroyd 6-constant fluid between concentric cylinders by the finite difference method. Hayat et al. 
[26] investigated the exact solution of a thin film flow of an Oldroyd 6 -constant fluid over a moving belt by the Homotopy 
perturbation. 

Investigation of nanofluid flow has received special focus in the past due to its relevance in numerous industrial applications. 
The researchers not only discovered unexpected thermal properties of nanofluids but also proposed mechanisms behind the 
enhanced thermal properties of nanofluids and thus identified unusual opportunities to develop them as next generation 
coolants for computers and safe coolants for nuclear reactors. A combination of nanofluid with biotechnological components 
can provide potential applications in agriculture, pharmaceuticals and biological sensors. Various types of nanomaterials 
including nanoparticles, nanowires, nanofibers, nanostructures, and nanomachines are used in biotechnological applications. 
The commercialization of nano-biotechnological products seems to have a potential future and within next a few years many 
new products of this nature are likely to be used. Nano and micro -fluidics is a new area which has potential for engineering 
applications, especially for the development of new biomedical devices and procedures [1], 
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The magnetic nanofluids possess both liquid and magnetic properties. These fluids have key importance in modulators, 
optical switches, optical gratings and tunable optical fiber filters. The magnetic nanoparticles are significant in medicine, 
construction of loud speakers, sink float separation, cancer therapy and tumor analysis. Nowadays, the sustainable energy 
generation is one of the most important issues across the globe. Perhaps the solar energy is also one of the best sources of 
renewable energy with minimal environment impact. Solar power is a direct way obtaining heat, water, and electricity from 
the nature. Researchers concluded that heat transfer and solar collection processes can be improved through the addition of 
nanoparticles in the fluids. Fluid heating and cooling are required in many industrial fields such as power manufacturing and 
transportation. Effective cooling techniques are needed for cooling any sort of high energy device. Common heat transfer 
fluids such as ethylene glycol, water and engine oil have limited heat transfer capabilities due to their low heat transfer 
properties and thus cannot meet the modem cooling requirements [1], 

The terminology of nanofluid was first used by Choi [22] when he experimentally discovered a significant improvement in 
the heat transfer performance through the addition of small amount nanometer sized particles in the base fluids. This addition 
also causes scattering of the incident radiation which allows a deeper absorption within the fluid. Recently, Trieb and Nitsch 
[3] proposed the idea of solar thermal collectors by using nanofluids to directly absorb the solar radiation. 

The phenomenon of thermal conductivity enhancement by dispersing nanoparticles in the liquid was also observed by 
Masuda et al. [7]. Buongiorno [9] recognized that the two main effects namely the Brownian motion and thermophoretic 
diffusion of nanoparticles contribute to the massive increase in the thermal conductivity of the liquids. He also proposed the 
modifications in the transport equations due to these effects. 

Vajravelu et al. [12] investigated heat and mass transfer properties of three -layer fluid flow in which nano-fluid layer is 
squeezed between two clear viscous fluid. Farooq et al. [27] studied heat and mass transfer of two-layer flows of third-grade 
nano-fluids in a vertical channel. Abou-zeid et al. [16] obtained numerical solutions and Global error estimation of natural 
convection effects on gliding motion of bacteria on a power-law nano-slime through a non-Darcy porous medium. El-Dabe et 
al. [17] investigated magneto-hydrodynamic non-Newtonian nano-fluid flow over a stretching sheet through a non-Darcy 
porous medium with radiation and chemical reaction. 

Magneto -hydrodynamic (MHD) is concerning the mathematical and physical scaffold that introduces magnetic -dynamics in 
electrically conducting fluids (e.g. in plasmas and liquid metals ). The applications of Magneto -hydrodynamic 
incompressible viscous flow in science and engineering involving heat and mass transfer under the influence of chemical 
reaction is great importance to many areas of science and engineering. This frequently occurs in agriculture, engineering, 
plasma studies and petroleum industries [15], 

Flow through porous media plays an important role in countless practical applications such as ground water flows, enhanced 
oil recovery processes, contamination of soils by hazardous wastes, pollution movement [20], An understanding of the 
dynamics of fluids in a porous medium is of principal interest because these flows are quite prevalent in nature. Such flows 
have attracted the attention of number of scholars due to their applications in many branches of science and technology, viz. 
in the fields of agricultural engineering to study the underground water resources; seepage of water in river beds; in 
petroleum technology to study the movement of natural gas, oil, and water through the oil reservoirs. In the widely used 
continuum approach to transport processes in a porous media, the differential equation governing the macroscopic fluid 
motion is based on the experimentally established Darcy law [8], which accounts for the drag exerted on the fluid by porous 
medium [29]. 

The objective of this work is to investigate the numerical solution by using Explicit Finite Difference method [6] for the 
system of non-linear differential equations which arises from magneto -hydrodynamic flows of an Oldroyd 6-constant Nano- 
fluid with a magnetic field through a non-Darcy porous medium with heat and mass transfer. We obtained the distributions of 
the velocity, the temperature, the concentration, and the Nanoparticles. Numerical results are found for different values of 
various non-dimensional parameters in the case of pure Coutte flow ( U 0 = 1, ). The effects of modified pressure 

gradients on those formulas were discussed in the case of pure Poiseuille flow and the generalized Couette flow. The results 
are shown graphically and discussed in detail. Also, the global error estimation for the error propagation is obtained by 
Zadunaisky technique [19], 

III. Fluid Model 

For an Oldroyd 6-constant fluid, the Cauchy stress tensor t is given by [30] 
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t=-PI+S *, (1) 

s *+h D ~^ + X i (S*A, + A,S*) + ^ (tr S*)A, = g [a ± + X 2 ^ + A 4 A\]. (2) 


Where, 


A 1 = L + L T , L — grad V . (3) 

—PI is the indeterminate part of the stress due to the constraint of incompressibility, S is the extra stress tensor, A 1 is the 

first Rivlin-Ericksen tensor, L is the velocity gradient, V is the velocity vector, and g , A ; (i = 1, 2, ,5) are the six material 

constants. 


The contra variant convected derivative ^ satisfying the principle of material frame in difference in terms of the material ^ 
is defined by the following equation 


DS* 

Dt 



s*l t . 


(4) 


in which the superscript T is used for the transpose operation. 


It should be noted that when = 0 (i = 1, 2, 3, 4, 5) the model reduces to the classical linearly viscous model of Newtonian 
fluid. For — 0 (i = 3, 4, 5) a 3-constant model of an Oldroyd-B fluid is described, while for A t = 0 (i = 2, 3, 4, 5) a 
Maxwell model is formulated. For A, = 0 (i = 1, 3, 4, 5) the model describes a second-grade fluid. 


IV. Mathematical Formulation 

We consider magnetohydrodynamic flows of an Oldroyd 6-constant Nano-fluid through a non-Darcy porous medium 
between two parallel plates (see FIG. 1). We take in our consideration the presence of heat and mass transfer with chemical 
reaction. 



(Oldroyed 6-Constant Nano- fluid) 

Figure 1: Magnetohydrodynamic flows of an Oldroyd 6-constant Nano-fluid through a non- 

Darcy porous medium between two parallel plates 

V. Basic Equations 

The basic equations governing the flow of an incompressible fluid are the following equations [10, 18], 

The continuity equation 

V.U = 0 , (5) 

The momentum equation 

Pf f = V - T + Ve (y A h) - PfC'V \V\ -^V + F , (6) 
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The energy equation 


The concentration equation 

dC 


(pc p ) / g] = k V 2 T + O + (pc p ) p [d b (VT. V0) + | (VT) 2 ] + 


Pf D B K T 

C, ’ 


dt 


= d b v 2 c + — v 2 r - a(c - c 2 y 


The nanoparticles concentration equation 


— = d b v 2 0 + — v 2 r 

dt u ^ T 2 


(7) 


( 8 ) 


(9) 


Where, ^ + (V. V)| denotes the material time derivative, J = erQ? + /r e F a//) , and E_— 0. 

For a simple one -dimensional shearing flow considered in this work, the extra stress tensor and the velocity are of the form 


' °xx °xy °XZ 


S(y, t) — Syx Syy S yz , l/(y ; t) — 


t $ZX S z y S zz J 


fu(y, t)N 
0 

V 0 > 


(10) 


Also, the constant magnetic field, the temperature, the concentration, and the nanoparticles are of the form 
H_ — (0, H 0 , 0) , T = T(y,t) , C = C(y,t), $ = $(y,t) . 


( 11 ) 


For steady motion, it is clear that with the simplified dependence (10-11), the continuity equation (5) is satisfied identically, 
and equations (1-3) and (6) become, 

£ = a ^ H o u ~ PC*u 2 ~^,u + p f ga T (T -T 2 ) + p f ga c (C - C 2 ) , (12) 


dP dSyy 
dy dy 

dP dS Z y 

dz dy 


S„ + (A, - 2A,)S„ 0 = M (h - 2A 2 ) g) Z 


S x y A-^Syy 


S’ + 


du 

dy 


S zy — = 0 

zy dy 


S yy + ^3 S xy = ^4 (^) 


S + —S 

*zy -f 2 -Jzx dy 


= 0 


S zz = o 


by addition of two equations (15) and (18) we have, 

+ Syy = 2 M U 4 - A 2 ) (gf - 2(A 3 - AdS xy g 
from equations (17) and (19) we get 


(13) 

(14) 

(15) 

(16) 

(17) 

(18) 

(19) 

( 20 ) 

( 21 ) 
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9 =9 =0 

' ’ z x - J zy 

Substituting (18) - (22) in (14) and (16) we have 



where, 

a l — ^1^-4 — C^-3 + 2-5 ) C-^-4 — ^ 2 ) 
a 2 = 2-i2-3 — (2.3 + A 5 )(/l3 — Ai) 


( 22 ) 

(23) 

(24) 

(25) 

(26) 


Defining the modified pressure 


P = P - S. 


yy 


and using (12) - (14) and (24) we arrive at 


dP 

dx 



— <7 nl Hq u — pc*u 2 


^ u + p f ga(T - T 2 ) + p f ga{C - C 2 ) 


dP 

dy 



(27) 


(28) 

(29) 


The dissipation function O can be written as follows 


dv i 

(D = T 

l J dx; 


By using the definitions (10) we have, 

du 


0 = s 


•xy 


dy 


Then, equations (7 - 9) can be written as follows: 

K fl + P _fEBKT^ + ( ) \ D + + 

dy 2 C s dy 2 P'p 15 dy dy 7* \dyj \ 


" dy 2 7i dy 2 v 


D B —j + — ~ 0 • 

" dy 2 r 2 dy 2 

The boundary conditions of the problem are given by: 

u = 0 , T = T X ,C = C 1 ,0 = 0 1 at y = 0 , 
u — u 0 , T — T 2 , C — C 2 , (p — 02 y = H . 


du fdu\ 

^ ^ “iWJ 


fdu\ 

i+ “ 2 W) 


(30) 


(31) 


— = 0 

(32) 

dy 


(33) 


(34) 


(35) 


We shall now write the field equations (28), (32 - 34) and the boundary conditions (35) in terms of a set of dimensionless 
variables and, for this purpose, we choose H and U as the characteristic length and velocity and introduce the following 
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dimensionless variables of both: 


y = - 

- / u 


u — - 
u 




CCi 


a 2 


( H /u ) 


a 7 — 7 ' X ~ h ’ ^ ~ 


C %) 


> K ' = 777 


0 = 


t-t 2 

Ti-T 2 


s = 


C—C 2 

c,-c 2 


0 


0-02 


(36) 


Thus, the system of our non-linear differential equations (28), (32 - 34) and the boundary conditions (35) can be rewritten in 
a non-dimensional forms: 


d 

dy 


du _ /duy 
dy + ai \dy) 


_ /duV 

1+ “ 2 W) 


dP n -2 

-- F s u z 

dx s 


( M + i) 


— ) H- G r Q H- B r S — 0 


1 d 2 0 , d 2 S , 

H - + iVj 

P r cLy 2 a dy 2 


*(f)(S) 


/de 


+ A H^) + 


du , _ (duy 
dy + ^Hdy/ 


1+ “ 2 (sf) 


= o , 


££!£ + Ld £!|_ i?c5m = 0 

L e dy 2 d dy 2 c 


d 2 0 , d 2 9 - 

Ah, — y + N t — - =0 

0 dy 2 c dy 2 


(37) 

(38) 

(39) 

(40) 


with the dimensionless boundary conditions: 

u — 0,0 = 1 ,5=1, 0 = 1 at y — 0 , 
u — U 0 , 6 — 0 ,5 = 0, 0 = 0 at y — 1 . 
Where, the dimensionless parameters are defined by 


(41) 


D„ = 


F, = 


pc*UH z 


, M — 


ap 2 H^H 2 _ga T H(T 1 -T 2 ) 


B r = 


g ac H (Ci-C 2 ) 


p _ M _ £>b/c7-(Ci-c 2 ) ^ _ d b{^i~4 > 2) _ Dt(T 1 -t 2 ) 

y /- ' / d f, /-rp rp -V / b t t 


k(Ti-T 2 ) 


9s c p v 0\ -T 2 ) 


Ti V 


, _ , _ Dt *T (n-T 2 ) _ 

e D B ' d T lV (Ci-C 2 ) ' c v 


(42) 


For convenience, we shall drop the bars that identify the dimensionless quantities. The system of non-linear ordinary 
differential equations (37) - (40) with the boundary conditions (41), will be solved numerically by using the explicit finite- 
difference method. And, we computed the global error for the solutions of the problem. 

VI. Numerical solution 

The equations (37) - (40) can be written after applied explicit finite difference schemes [6] as: 


( 




tip— 1] — 2ti[t]+u[i + l]\ l^ +Fs( “ W)2+ ( M+ ^) u[t] ~ Gr 9li] - Br Sp]].[l+«2(“ [ ‘ ] h lL l! ) 


= 0 


^e[i-i]-2e[i]+e[i+i]^ 


1 + 3 ai (!Mizll) _ a 2 (M=ll) +aia 2 («H) j 

^ 5p-l]-25[i]+5[i+l] ^ + p N ^ ^ flp]-flp-l] ^ ^ 0[l]-0[i-l] j + p N ^ ^ flpl-flp-1] ^ 


(43) 


I + Pr N d 


1 +gl (HlWl) 


1+a2 (HWz^z 1!) 


^ u[i]-u[i-l] j 2 _ q 


+ 

(44) 
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(Sli- 


l]-2S[i]+S[i + l] 


h 2 


) + L e L d (— 


l]-2Q[i]+6[i+l] 


h z 


) - L e R c (S[i]) m = 0 , 


0U-i]-20[i]+0[i+i] j + N t j ^ eu-i]-2e[i]+e[t+i] '\ _ 


h 2 


h 2 


) = °- 


(45) 

(46) 


Where the index i refers to y and the A y — h — 0.04 . According to the boundary conditions (41) we can solved equations 
(43) - (46) numerically, then a Newtonian iteration method continues until either of goals specified by accuracy goal or 
precision goal is achieved. 


VII. Estimation of the global error 

We used Zadunaisky technique [19] for calculating the global error, which can be explained in the following steps: 

1. We fined the interpolating polynomial of u, u, 9,9, S, S ', (f), (/) from the values of them which came from the 

explicit finite -difference method. Then we named it P[{x) , (i = 1,2,3, ,8) , and we found the interpolating 

functions of u ', O ' , S ', and <p , and we named them as: 

u" = q^x), 9" = q 2 {x), S" = q 3 (x), cp" = q 4 (x) , 

^+F s (p 1 (X)) 2 +(m+A-) P 1 (x)-G r P 3 (x)-B r p 5 (x)].[l+a 2 (p 2 (x)) 2 
[l+3a 1 (P 2 (x)) 2 -a 2 (p 2 (x)) 2 +a 1 a: 2 (P 2 (x)) 4 ] 


<7l00 



<72 O) = 


1 

(l — Nd Le 


-P r N d L e R c {P s {x)) m 


P r N b P 4 (x) P B (x) - P r N t (P 4 (x)) 2 - P r 


l+ g i(P 2 Qr)) 2 

X+a 2 (P 2 (x )) 2 


(PzW ) 2 


<?3 W = -L e Ld QzM + L e R c (P 5 (x)) m 


q 4 (x) = - Nt / Nh q 2 (x)- (47) 

2. We calculate the defect functions D[(x) , (i = 1,2,3, ,8) , which can be written as follows: 

£>i(x) = P[ (x) - P 2 {x) = 0 , D 2 (x) = P" (x) - q 4 (x) , 

D 3 W = P 3 W - P 4 (x) = 0 , £> 4 (x) = P 3 "(x) - q 2 {x ) , 

D 5 W = -P 6 (x) = 0 , £> 6 (x) = Ps'(x) - q 3 (x) , 

£>7 to = £7 to - P&(x) = 0 , £>g Or) = P"(x) - q 4 (x) . (48) 

3. We add the defect functions Dj(x) , (i = 1,2,3, ,8) to the original problem. 

4. We solved the pseudo-problem (new problem) by the same method which is used for solving the main problem, and we 
will have the new solutions. 

5. We calculate an estimation of the global error from the formulas, 

e„ = /„ - f_{Xn) = fn~ P(x n ) , (n = 0,1,2, ,25) (49) 

In this relation, f n is the approximate solutions of the new problem (the pseudo -problem) at the point x n , and f(x n ) is the 
exact solutions of pseudo-problem at x n . 

The values of global error for the solutions of the problem which solved by the explicit finite difference method are shown in 
table (1). The error in table (1) based on using 26 points to find interpolating polynomials P;(x) of degree 25. In order to 
achieve the above task, we used the Mathematica package. 
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Table 1 

Different values of the dimensionless physical quantities 


The global error by using finite difference technique 

y 

u =Y1 

error (eln) 

e=Y3 

error (e3n) 

S=Y5 

error (e Sn ) 

<t»=Y7 

error (e7n) 

0 

0.00E+00 

0.00E+00 

1.00E+00 

0.00E+00 

1.00E+00 

0.00E+00 

1.00E+00 

0.00E+00 

0.08 

0.084677 

2.45938E-06 

1.0654 

1.06548E-05 

0.877711 

4.73414E-06 

0.832759 

9.4549 IE-06 

0.16 

0.167804 

4.11468E-06 

1.10363 

1.91716E-05 

0.765806 

8.26993E-06 

0.681824 

1.69842E-05 

0.24 

0.249488 

5.13257E-06 

1.11623 

2.58221E-05 

0.662794 

1.08326E-05 

0.546264 

2.28215E-05 

0.32 

0.329882 

5.63641E-06 

1.10429 

3.07583E-05 

0.567513 

1.25695E-05 

0.425427 

2.71033E-05 

0.4 

0.40918 

5.70977E-06 

1.06848 

3.40436E-05 

0.479077 

1.35718E-05 

0.318913 

2.98941E-05 

0.48 

0.487619 

5.39798E-06 

1.00906 

3.56739E-05 

0.39684 

1.38901E-05 

0.226562 

3.12032E-05 

0.56 

0.565479 

4.70768E-06 

0.925913 

3.55924E-05 

0.320364 

1.3545E-05 

0.148452 

3.09962E-05 

0.64 

0.64308 

3.60483E-06 

0.818507 

3.37017E-05 

0.249402 

1.25357E-05 

0.0848959 

2.92048E-05 

0.72 

0.720789 

2.019E-06 

0.685922 

2.98767E-05 

0.183879 

1.0847E-05 

0.0364469 

2.57363E-05 

0.8 

0.799017 

9.48686E-08 

0.526812 

2.39868E-05 

0.123881 

8.45801E-06 

0.00391268 

2.04927E-05 

0.88 

0.878223 

2.40736E-06 

0.339383 

1.59365E-05 

0.0696485 

5.37175E-06 

-0.0116299 

1.34196E-05 

0.96 

0.958919 

2.76915E-06 

0.121352 

5.76248E-06 

0.0215698 
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VIII. Numerical results and discussion 

In this paper we generalized the problem of MHD flows to include the non-Newtonian Nano-fluid obeying Oldroyd 6- 
constant model through a porous medium of non-Darcy type with heat and mass transfer. The system of non-linear ordinary 
differential equations (37) - (40) with the boundary conditions (41) was solved numerically by using an Explicit finite 
Difference method. The functions u ,6, S, and <p are obtained and illustrated graphically as shown in figures (a-1 - a-5), 
(b-1 - b-13), (c-1 - c-13), and (d-1 - d-13) for different values of the parameters of the problem in the case of pure Couette 
flow (U 0 = 1, ^=o) , at the values of non-Newtonian parameters (oq = 8, oq = 2). And, figures (e-1 - e-4) show the 
effects of modified pressure gradients on those formula in the case of pure Poiseuille flow ( U 0 — 0), at the values of non- 
Newtonian parameters (oq = 8, a 2 — 2). And, figures (f-1 - f-4) show the effects of modified pressure gradients on those 
formula in the case of generalized Couette flow (U 0 = 1), at the values of non-Newtonian parameters (oq = 8, a 2 — 2). 



FIG. A-1: Profiles of the velocity u(y) for a pure Couette flow with various values of Br for a system have 
the particulars M=0.1, Da=l, Fs=0.5, Gr=0.1, Pr=0.7, Nd=2, Nb=0.5, Nt=0.3, Le=2, Ld=0.1, Rc=0.5, m=2. 
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FIG. A-2: Profiles of the velocity u(y) for a pure Couette flow with various values of Gr for a system have 
the particulars M=0.1, Da=l, Fs=0.5, Br=0.1, Pr=0.7, Nd=2, Nb=0.5, Nt=0.3, Le=2, Ld=0.1, Rc=0.5, m=2. 



FIG. A-3: Profiles of the velocity u(y) for a pure Couette flow with various values of Da for a system have 
the particulars M=0.1, Fs=0.5, Nt=0.3, Gr=0.1, Br= 0.1, Pr=0.7, Nd=2, Nb=0.5, Nt=0.3, Le=2, Ld=0.1, 

Rc=0.5, m=2. 
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FIG. A-4: Profiles of the velocity u(y) for a pure Couette flow with various values of Fs for a system have 
the particulars M=0.1, Da=l, Nt=0.3, Gr=0.1, Br=0.1, Pr=0.7, Nd=2, Nb=0.5, Nt=0.3, Le=2, Ld=0.1, Rc=0.5, 

m=2. 
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FIG. A-5: Profiles of the velocity u(y) for a pure Couette flow with various values of M for a system have 
the particulars Da=l, Fs=0.5, Nt=0.3, Gr=0.1, Br= 0.1, Pr=0.7, Nd=2, Nb=0.5, Nt=0.3, Le=2, Ld=0.1, 

Rc=0.5, m=2. 


8.1 Couette flow 

Considering the Couette flow of an Oldroyd 6-constant fluid between two infinite parallel plates, for which the lower plate is 
fixed and the top plate moves with the velocity (U 0 — 1), it is assumed that the flow is driven only by the motion of the top 
plate. The modified pressure gradient in the flow direction vanishes instant (^|=o). Figures (a-1 - a-5), (b-1 - b-13), (c-1 - c- 
13), and (d-1 - d-13) will show the effect of the problem parameters on the problem solutions in this case. 

Figures (a-1 - a-5) show the distributions of the velocity profile u at different values of some parameters of the problem. It is 
clear that, the velocity increases by increasing each of local temperature Grashof number G r , local nanoparticle Grashof B r , 
and Darcy number D a . But, the velocity decreases by increasing each of Forchheimer number F s , magnetic field parameter M. 

Figures (b-1 - b-13) describe the distributions of the temperature profile 0 at different values of some parameters of the 
problem. It is noted that, as the local nanoparticle Grashof B r , local temperature Grashof number G r , and Darcy number D a 
increase, the temperature increases in the region 0 < y < 0.58 and, it returns decrease to 0.58 < y < 1. The temperature 
increases by increasing each of Sort number L d , Lewis number L e , Brownian motion parameter N b , Dufour number N d , 
Thermophoresis parameter N t , Prandtl number P r , and Chemical reaction parameter R c . And, as Forchheimer number F s 
increases, the temperature decreases in the region 0 < y < 0.7 and, it returns increase to 0.7 < y < 1. Also, as magnetic 
field parameter M increases, the temperature decreases in the region 0 < y < 0.5 and, it returns increase to 0.5 < y < 1. 
Finally, the temperature decreases by increasing Chemical reaction order m. 



FIG. B-1: Profiles of the temperature 0(y) for a pure Couette flow with various values of Br for a system 
have the particulars M=0.1, Da=l, Fs=0.5, Gr=0.1, Pr=0.7, Nd=2, Nb=0.5, Nt=0.3, Le=2, Ld=0.1, Rc=0.5, 

m=2. 
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FIG. B-2: Profiles of the temperature 0(y) for a pure Couette flow with various values of Gr for a system 
have the particulars M=0.1, Da=l, Fs=0.5, Br=0.1, Pr=0.7, Nd=2, Nb=0.5, Nt=0.3, Le=2, Ld=0.1, Re=0.5, 

m=2. 



FIG. B-3: Profiles of the temperature 0(y) for a pure Couette flow with various values Ld for a system 
have the particulars M=0.1, Da=l, Fs=0.5, Gr=0.1, Br=0.1, Pr=0.7, Nd=2, Nb=0.5, Nt=0.3, Le=2, Rc=0.5, 

m=2. 



FIG. B-4: Profiles of the temperature 0(y) for a pure Couette flow with various values of Le for a system 
have the particulars M=0.1, Da=l, Fs=0.5, Gr=0.1, Br=0.1, Pr=0.7, Nd=2, Nb=0.5, Nt=0.3, Ld=0.1, Rc=0.5, 

m=2. 
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FIG. b-5: Profiles of the temperature 0(y) for a pure Couette flow with various values of Nb for a system 
have the particulars M=0.1, Da=l, Fs=0.5, Gr=0.1, Br=0.1, Pr=0.7, Nd=2, Le=2, Nt=0.3, Ld=0.1, Re=0.5, 

m=2. 



FIG. b-6: Profiles of the temperature 0(y) for a pure Couette flow with various values of Nd for a system 
have the particulars M=0.1, Da=l, Fs=0.5, Gr=0.1, Br=0.1, Pr=0.7, Nb=0.5, Le=2, Nt=0.3, Ld=0.1, Rc=0.5, 

m=2. 



FIG. b-7: Profiles of the temperature 0(y) for a pure Couette flow with various values of Nt for a system 
have the particulars M=0.1, Da=l, Fs=0.5, Gr=0.1, Br=0.1, Pr=0.7, Nb=0.5, Le=2, Nd=2, Ld=0.1, Rc=0.5, 

m=2. 
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FIG. b-8: Profiles of the temperature 0(y) for a pure Couette flow with various values of Da for a system 
have the particulars M=0.1, Fs=0.5=l, Nt=0.3, Gr=0.1, Br=0.1, Pr=0.7, Nd=2, Nb=0.5, Nt=0.3, Le=2, 

Ld=0.1, Rc=0.5, m=2. 
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FIG. b-9: Profiles of the temperature 0(y) for a pure Couette flow with various values of Fs for a system 
have the particulars M=0.1, Da=l, Nt=0.3, Gr=0.1, Br=0.1, Pr=0.7, Nd=2, Nb=0.5, Nt=0.3, Le=2, Ld=0.1, 

Rc=0.5, m=2. 



FIG. b-10: Profiles of the temperature 0(y) for a pure Couette flow with various values of M for a system 
have the particulars Da=l, Fs=0.5=l, Nt=0.3, Gr=0.1, Br=0.1, Pr=0.7, Nd=2, Nb=0.5, Nt=0.3, Le=2, Ld=0.1, 

Rc=0.5, m=2. 
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FIG. b-11: Profiles of the temperature 0(y) for a pure Couette flow with various values of Pr for a system 
have the particulars Da=l, Fs=0.5=l, Nt=0.3, Gr=0.1, Br=0.1, M=0.1, Nd=2, Nb=0.5, Nt=0.3, Le=2, Ld=0.1, 

Rc=0.5, m=2. 



FIG. b-12: Profiles of the temperature 0(y) for a pure Couette flow with various values of Rc for a system 
have the particulars Pr=0.7, Da=l, Fs=0.5=l, Nt=0.3, Gr=0.1, Br=0.1, M=0.1, Nd=2, Nb=0.5, Nt=0.3, Le=2, 

Ld=0.1, m=2. 



FIG. b-13: Profiles of the temperature 0(y) for a pure Couette flow with various values of m for a system 
have the particulars Pr=0.7, Da=l, Fs=0.5=l, Nt=0.3, Gr=0.1, Br=0.1, M=0.1, Nd=2, Nb=0.5, Nt=0.3, Le=2, 

Ld=0.1, Rc=0.5. 
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Figures (c-1 - c-13) illustrate the distributions of the concentration profile S at different values of some parameters of the 
problem. It is seen that, the concentration decreases by increasing each of local nanoparticle Grashof B r , local temperature 
Grashof number G r , Sort number L d , Lewis number L e , Brownian motion parameter N b , Dufour number N d , Thermophoresis 
parameter N t , Prandtl number P r , and Chemical reaction parameter R c . And, as Darcy number D a increases, the concentration 
decreases in the region 0 < y < 0.6 and, it returns increase to 0.6 < y < 1. But, as Forchheimer number F s increases, the 
concentration increases in the region 0 < y < 0.7 and, it returns decrease to 0.7 < y < 1. The effect of magnetic field 
parameter M on the concentration disappears in the region 0 < y < 0.6, the concentration returns decrease to 0.6 < y < 1 at 
large values of M. Finally, the concentration increases by increasing Chemical reaction order m. 



FIG. c-1: Profiles of the concentration S(y) for a pure Couette flow with various values of Br for a system have the 
particulars M=0.1, Da=l, Fs=0.5, Gr=0.1, Pr=0.7, Nd=2, Nb=0.5, Nt=0.3, Le=2, Ld=0.1, Rc=0.5, m=2. 



FIG. c-2: Profiles of the concentration S(y) for a pure Couette flow with various values of Gr for a system 
have the particulars M=0.1, Da=l, Fs=0.5, Br=0.1, Pr=0.7, Nd=2, Nb=0.5, Nt=0.3, Le=2, Ld=0.1, Rc=0.5, 

m=2. 



FIG. c-3: Profiles of the concentration S(y) for a pure Couette flow with various values of Ld for a system have the 
particulars M=0.1, Da=l, Fs=0.5, Gr=0.1, Br=0.1, Pr=0.7, Nd=2, Nb=0.5, Nt=0.3, Le=2, Rc=0.5, m=2. 
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FIG. c-4: Profiles of the Concentration S(y) for a pure Couette flow with various values of Le for a system have the 
particulars M=0.1, Da=l, Fs=0.5, Gr=0.1, Br=0.1, Pr=0.7, Nd=2, Nb=0.5, Nt=0.3, Ld=0.1, Rc=0.5, m=2. 



FIG. c-5: Profiles of the concentration S(y) for a pure Couette flow with various values of Nb for a system have the 
particulars M=0.1, Da=l, Fs=0.5, Gr=0.1, Br=0.1, Pr=0.7, Nd=2, Le=2, Nt=0.3, Ld=0.1, Rc=0.5, m=2. 



FIG. c-6: Profiles of the Concentration S(y) for a pure Couette flow with various values of Nd for a system have the 
particulars M=0.1, Da=l, Fs=0.5, Gi-0.1, B r=0.1, Pr=0.7, Nb=0.5, Le=2, Nt=0.3, Ld=0.1, Rc=0.5, m=2. 



FIG. c-7: Profiles of the Concentration S(y) for a pure Couette flow with various values of Nt for a system have the 
particulars M=0.1, Da=l, Fs=0.5, Gr=0.1, Br=0.1, Pr=0.7, Nb=0.5, Le=2, Nd=2, Ld=0.1, Rc=0.5, m=2. 
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FIG. c-8: Profiles of the concentration S(y) for a pure Couette flow with various values of Da for a system have the 
particulars M=0.1, Fs=l, Nt=0.3, Gr=0.1, Br=0.1, Pr=0.7, Nd=2, Nb=0.5, Nt=0.3, Le=2, Ld=0.1, Rc=0.5, m=2. 



FIG. c-9: Profiles of the concentration S(y) for a pure Couette flow with various values of Fs for a system have the 
particulars M=0.1, Da=l, Nt=0.3, Gr=0.1, Br=0.1, Pr=0.7, Nd=2, Nb=0.5, Nt=0.3, Le=2, Ld=0. 1, Rc=0.5, m=2. 



FIG. c-10: Profiles of the concentration S(y) for a pure Couette flow with various values of M for a system have the 
particulars Da=l, Fs=l, Nt=0.3, Gi-0.1, Br=0.1, Pr=0.7, Nd=2, Nb=0.5, Nt=0.3, Le=2, Ld=0.1, Rc=0.5, m=2. 
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FIG. c-11: Profiles of the concentration S(y) for a pure Couette flow with various values of Pr for a system have the 
particulars Da=l, Fs=l, Nt=0.3, Gr=0.1, Br=0.1, M=0.1, Nd=2, Nb=0.5, Nt=0.3, Le=2, Ld=0.1, Rc=0.5, m=2. 



FIG. c-12: Profiles of the concentration S(y) for a pure Couette flow with various values of Rc for a system have the 
particulars Pr=0.7, Da=l, Fs=l, Nt=0.3, Gr=0.1, Br=0.1, M=0.1, Nd=2, Nb=0.5, Nt=0.3, Le=2, Ld=0.1, m=2. 



FIG. c-13: Profiles of the concentration S(y) for a pure Couette flow with various values of m for a system have the 
particulars Pr=0.7, Da=l, Fs=l, Nt=0.3, Gr=0.1, Br=0.1, M=0.1, Nd=2, Nb=0.5, Nt=0.3, Le=2, Ld=0.1, Rc=0.5. 

Figures (d-1 - d-13) display the distributions of the nanoparticle profile <p at different values of some parameters of the 
problem. It is clear that, , the nanoparticle decreases by increasing each of local nanoparticle Grashof B r , local temperature 
Grashof number G r , Sort number L d , Lewis number L e , Dufour number N d , Thermophoresis parameter N t , Prandtl number P r , 
and Chemical reaction parameter R c . But, the nanoparticle increases by increasing Brownian motion parameter N b And, as 
Darcy number D a increases, the nanoparticle decreases in the region 0 < y < 0.57 and, it returns increase to 0.57 < y < 1. 
But, as Forchheimer number F s increases, the nanoparticle increases in the region 0 < y < 0.7 and, it returns decrease to 
0.7 < y < 1. The effect of magnetic field parameter M on the nanoparticle disappears in the region 0 < y < 0.5 and, the 
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nanoparticle returns decrease to 0.5 < y < 1 at large values of M. Finally, the nanoparticle increases by increasing Chemical 
reaction order m. 



FIG. d-1: Profiles of the nano-particle <f>(\ ) for a pure Couette flow with various values of Br for a system have the 
particulars M=0.1, Da=l, Fs=0.5, Gr=0.1, Pr=0.7, Nd=2, Nb=0.5, Nt=0.3, Le=2, Ld=0.1, Rc=0.5, m=2. 



FIG. d-2: Profiles of the nano-particle </)( y) for a pure Couette flow with various valuesof Gr for a system have the 
particulars M=0.1, Da=l, Fs=0.5, Bi-0.1, Pr=0.7, Nd=2, Nb=0.5, Nt=0.3, Le=2, Ld=0.1, Rc=0.5, m=2. 



FIG. d-3: Profiles of the nano-particle 0(y) for a pure Couette flow with various values of Ld for a system have the 
particulars M=0.1, Da=l, Fs=0.5, Gr=0.1, Br=0.1, Pr=0.7, Nd=2, Nb=0.5, Nt=0.3, Le=2, Rc=0.5, m=2. 
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FIG. d-4: Profiles of the nano-particle 0(v) for a pure Couette flow with various values the Le for a system have the 
particulars M=0.1, Da=l, Fs=0.5, Gr=0.1, Br=0.1, Pr=0.7, Nd=2, Nb=0.5, Nt=0.3, Ld=0.1, Rc=0.5, m=2. 



FIG. d-5: Profiles of nano-particle <f>( y) for a pure Couette flow with various values of Nb for a system have the 
particulars M=0.1, Da=l, Fs=0.5, Gr=0.1, Br=0.1, Pr=0.7, Nd=2, Le=2, Nt=0.3, Ld=0.1, Rc=0.5, m=2. 



FIG. d-6: Profiles of the nano-particle <f>(\) for a pure Couette flow with various values of Nd for a system have the 
particulars M=0.1, Da=l, Fs=0.5, Gr=0.1, Br=0.1, Pr=0.7, Nb=0.5, Le=2, Nt=0.3, Ld=0.1, Rc=0.5, m=2. 
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FIG. d-7: Profiles of the nano-particle for a pure Couette flow with various values of Nt for a system have the 
particulars M=0.1, Da=l, Fs=0.5, Gr=0.1, Br=0.1, Pr=0.7, Nb=0.5, Le=2, Nd=2, Ld=0.1, Rc=0.5, m=2. 
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FIG. d-8: Profiles of the nano-particle <f>(\ ) for a pure Couette flow with various values of Da for a system have the 
particulars M=0.1, Fs=0.5, Nt=0.3, Gr=0.1, Br=0.1, Pr=0.7, Nd=2, Nb=0.5, Nt=0.3, Le=2, Ld=0.1, Rc=0.5, m=2. 



FIG. d-9: Profiles of the nano-particle <p( y) for a pure Couette flow with various values of Fs for a system have the 
particulars M=0.1, Da=l, Nt=0.3, Gr=0.1, Br=0.1, Pr=0.7, Nd=2, Nb=0.5, Nt=0.3, Le=2, Ld=0.1, Rc=0.5, m=2. 
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FIG. d-10: Profiles of the nano-particle <p(y) for a pure Couette flow with various values of M for a system have the 
particulars Da=l, Fs=0.5, Nt=0.3, Gr=0.1, Br=0.1, Pr=0.7, Nd=2, Nb=0.5, Nt=0.3, Le=2, Ld=0.1, Rc=0.5, m=2. 



FIG. d-11: Profiles of the nano-particle <p(y) for a pure Couette flow with various values of Pr for a system have the 
particulars Da=l, Fs=0.5, Nt=0.3, Gr=0.1, Br=0.1, M=0.1, Nd=2, Nh=0.5, Nt=0.3, Le=2, Ld=0.1, Rc=0.5, m=2. 



FIG. d-12: Profiles of the nano-particle <p{y) for a pure Couette flow with various values of Rc for a system have the 
particulars Pr=0.7, l)a=l, I s=0.5, Nt=0.3, Gr=0.1, Br=0.1, M=0.1, Nd=2, Nb=0.5, Nt=0.3, Le=2, Ld=0.1, m=2. 
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FIG. d-13: Profiles of the nano-particle <f)(y) for a pure Couette flow with various values of m for a system have the 
particulars Pr=0.7, Da=l, Fs=0.5, Nt=0.3, Gr=0.1, Br=0.1, M=0.1, Nd=2, Nb=0.5, Nt=0.3, Le=2, Ld=0.1, Rc=0.5. 


8.2 Poiseuille flow 

d.P 

Figures (e-1 - e-4) describe the effect of pressure gradient ( — ) on the velocity, the temperature, the concentration, and the 

dP 

nanoparticle distributions. It is clear that, the velocity, the concentration, and the nanoparticle decrease with increasing ( — ), 

dP 

while the temperature increases with increasing ( — ). 



FIG. e-1: Profiles of the velocity u(y) for a pure Poiseuille flow with various values of dP/dx for a system have the 
particulars M=0.1, Da=l, Fs=0.5, Gr=0.1, Br=0.1, Pr=0.7, Nd=2, Nb=0.5, Nt=0.3, Le=2, Ld=0.1, Rc=0.5, m=2. 



FIG. e-2: Profiles of the temperature 0(y) for a pure Poiseuille flow with various values of dP/dx for a system have 
the particulars M=0.1, Da=l, Fs=0.5, Gr=0.1, Br=0.1, Pr=0.7, Nd=2, Nb=0.5, Nt=0.3, Le=2, Ld=0.1, Rc=0.5, m=2. 
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FIG. e-3: Profiles of the concentration S(y) for a pure Poiseuille flow with various values of dP/dx for a system have 
the particulars M=0.1, Da=l, Fs=0.5, Gr=0.1, Br=0.1, Pr=0.7, Nd=2, Nb=0.5, Nt=0.3, Le=2, Ld=0.1, Rc=0.5, m=2. 



FIG. e-4: Profiles of the nano-particle </>(>) for a pure Poiseuille flow with various values of dP/dx for a system have 
the particulars M=0.1, Da=l, Fs=0.5, Gr=0.1, Br=0.1, Pr=0.7, Nd=2, Nb=0.5, Nt=0.3, Le=2, Ld=0.1, Rc=0.5, m=2. 

8.3 Generalized Couette flow 

If the flow of an Oldroyd 6-constant fluid between two parallel plates is driven by both the motion of the top plate and a 

dP 

constant pressure gradient in the direction parallel to the plates, a generalized couette flow is formed for ( — yO, U 0 =l)- 
Numerical results for such a generalized couette flow are illustrated in figures (f-1 - f-4). 



FIG. f-1: Profiles of the velocity u(y) for a generalized Couette flow with various values of dP/dx for a system have 
the particulars M=0.1, 1)a=l, Fs=0.5, Gr=0.1, Br=0.1, Pr=0.7, Nd=2, Nb=0.5, Nt=0.3, Le=2, 1x1=0. 1, Rc=0.5, m=2. 
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FIG. f-2: Profiles of the temperature 0(y) for a generalized Couette flow with various values of dP/dx for a system 
have the particulars M=0.1, Da=l, Fs=0.5, Gr=0.1, Br=0.1, Pr=0.7, Nd=2, Nb=0.5, Nt=0.3, Le=2, Ld=0.1, Rc=0.5, m=2. 



FIG. f-3: Profiles of the concentration S(y) for a generalized Couette flow with various values of dP/dx for a system 
have the particulars M=0.1, Da=l, Fs=0.5, Gr=0.1, Br=0.1, Pr=0.7, Nd=2, Nb=0.5, Nt=0.3, Le=2, Ld=0.1, Rc=0.5, m=2. 



FIG. f-4: Profiles of the nano-particle </>(y) for a generalized Couette flow with various values of dP/dx for a system 
have the particulars M=0.1, Da=l, Fs=0.5, Gr=0.1, Br=0.1, Pr=0.7, Nd=2, Nb=0.5, Nt=0.3, Le=2, Ld=0.1, Rc=0.5, m=2. 

Figures (f-1 - f-4) display the variation of the velocity, the temperature, the concentration, and the nanoparticle distributions 

cLP 

for several values of the pressure gradient for ( — ). It is observed that, the velocity, the concentration, and the nanoparticle 

cLP cLP 

decrease with increasing ( — ), while the temperature increases with increasing ( —). 


International Journal of Engineering Research & Science (IJOER) 


ISSN, [2395-6992] 


[Vol-3, Issue-7, July- 2017] 


IX. Conclusion 

In this work, we have studied magnetohydrodynamic flows of an Oldroyd 6-constant nano-fluid through a non-Darcy porous 
medium with heat and mass transfer. The governing boundary value problem was solved numerically by an Explicit Finite - 
Difference method. We concentrated our work on obtaining the velocity, the temperature, the concentration, and the 
nanoparticle distributions which are illustrated graphically at different values of the parameters of the problem in three cases, 
pure Couette flow, pure Poiseuille flow, and generalized Couette flow. Global error estimation is also obtained using 
Zadunaisky technique. We used 26 points to find the interpolating polynomial of degree 25 in interval [0,1] and the results 
are shown in table (1). We notice that, the error in table (1) is good enough to justify the use of resulting numerical values. 
The main findings from the current study can be summarized as follows: 

1- By increasing B r , G r , and D a , the velocity u increases, whereas it decreases by increasing each of F s , and M in the case 
of Couette flow at (a l =8 ,a 2 = 2). And, the other parameters doesn’t have any effect on the velocity. 

2- By increasing B r , G r , and D a , the temperature 0 increases at the first part (0 < y < 0.58) while, it decreases at the 
second part (0.58 < y < 1). And, by increasing L d , L e , N b , N d , N t , P r , and R c the temperature increases, while by 
increasing m, the temperature decreases. By increasing F s , the temperature decreases at the first part (0 < y < 0.7) 
while, it increases at the second part (0.7 < y < 1). Also, By increasing M, the temperature decreases at the first part 
(0 < y < 0.5) while, it increases at the second part ( 0.5 < y < 1). All these results are given in the case of Couette 
flow at ( a 1 =8 ,a 2 = 2). 

3- By increasing Br, Gr, Ld, Le, Nb, Nd, Nt, Pr, and Rc, the concentration S decreases, while by increasing m, the 
concentration increases. By increasing Da, the concentration decreases at the first part (0 < y < 0.6) while, it 
increases for (0.6 < y < 1). But, by increasing Fs, the concentration increases in the first part (0 < y < 0.7) while, it 
decreases at the second part (0.7 < y < 1). And, by increasing M, the concentration doesn’t have any change at the 
first part (0 < y < 0.6) while, it decreases at the second part (0.6 < y < 1) at large values of M. All these results are 
given in the case of Couette flow at (aq = 8, a 2 — 2). 

4- By increasing B r , G r , L d , L e , N d , N t , P r and R c , the nanoparticle (f) decreases, while by increasing N b , and m, the 
nanoparticle increases. By increasing Da, the nanoparticle decreases at the first part (0 < y < 0.57) while, it increases 
for (0.57 < y < 1). But, by increasing Fs, the nanoparticle increases in the first part (0 < y < 0.7) while, it decreases 
at the second part (0.7 < y < 1). And, by increasing M, the nanoparticle doesn’t have any change at the first part 
(0 < y < 0.5) while, it decreases at the second part (0.5 < y < 1) at large values of M. All these results are given in 
the case of Couette flow at (aq — 8, a 2 =2). 

dp 

5- By increasing — , the velocity, the concentration, and the nanoparticle decrease, while the temperature increases in 
the case of Poiseuille flow at (aq = 8, a 2 = 2). 

dP 

6- By increasing — , the velocity, the concentration, and the nanoparticle decrease, while the temperature increases in 
the case of generalized Couette flow at ( a 1 = 8 ,a 2 =2). 
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